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Abstract:  In  a  previous  paper  we  developed  a  matrix  theory  that 
applies  to  any  cylindrically  symmetric  fiber  surrounded  by  Bragg 
cladding.  Using  this  formalism,  along  with  Finite  Difference  Time  Do¬ 
main  (FDTD)  simulations,  we  study  the  waveguide  dispersion  for  the 
m  =  1  mode  in  an  air-core  Bragg  fiber  and  showed  it  is  possible  to 
achieve  very  large  negative  dispersion  values  (~  -20,000  ps/(nm.km)) 
with  significantly  reduced  absorption  loss  and  non-linear  effects. 
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1  Introduction 

The  idea  of  confining  light  in  a  Bragg  fiber  was  first  proposed  by  Yeh  et.  al.  in  [1].  Unlike 
conventional  fibers,  which  guide  light  through  total  internal  reflection,  Bragg  fibers 
achieve  light  confinement  through  Bragg  reflection.  Since  Bragg  fibers  and  conventional 
optical  fibers  utilize  different  guiding  mechanisms,  it  is  not  surprising  that  Bragg  fibers 
offer  many  interesting  features  that  are  difficult  to  achieve  in  conventional  fibers.  One 
such  example  is  the  possibility  of  guiding  light  in  air,  which  has  attracted  much  recent 
interest  [2]- [10].  In  a  Bragg  fiber  with  large  cladding  index  contrast,  light  can  be  confined 
mostly  within  the  center  air  core,  which  can  lead  to  lower  absorption  loss  and  reduce  the 
threshold  for  nonlinear  effects  [5] .  Another  area  of  particular  interest  is  the  dispersion 
characteristics  of  the  air-core  Bragg  fiber.  Johnson  et  al.  [9]  have  shown  that  by  making 
the  air  core  sufficiently  large,  arbitrarily  small  group  velocity  dispersions  can  be  achieved 
for  the  TEoi  mode  (D  ~  12ps/(nm.km)  has  been  demonstrated),  as  dispersion  D  in 
this  case  scales  as  l/r%ore,  where  rcore  is  the  core  radius.  In  this  paper  we  focus  on 
a  different  application:  dispersion  compensation,  for  which  a  large,  negative  dispersion 
value  is  desired.  Apparently  a  different  approach  should  be  used  in  our  case,  since  large 
dispersion  values  are  generally  associated  with  small  core  radii. 

For  years  researchers  have  proposed  various  schemes  to  compensate  chromatic  dis¬ 
persion  using  DCFs  [11]- [15] ,  with  dispersion  values  D  typically  ranging  up  to  -1,000 
ps/(ns.km).  Since  none  of  the  above-mentioned  schemes  guide  light  in  air,  they  in¬ 
evitably  suffer  from  absorption  loss  and  nonlinearity  effects  such  as  self-phase  modula¬ 
tion  and  four- wave  mixing  [16].  The  degree  to  which  these  undesirable  effects  undermine 
the  performance  of  the  DCFs  depends  on  each  particular  design.  In  the  following  sec¬ 
tions  we  will  show  that  it  is  theorectically  possible  to  design  a  Bragg  fiber  that  not  only 
guides  light  in  air,  but  also  produces  large  negative  dispersions  (~  -20,000  ps/(nm.km)). 

2  Theoretical  analysis 

In  Ref.  [10]  we  developed  a  general  matrix  theory  for  cylindrically  symmetric  waveguides 
with  Bragg  claddings.  In  this  theory,  an  arbitrary  number  of  inner  dielectric  layers, 
called  the  “inner  region”  (IR)1,  are  treated  exactly,  whereas  the  outside  cladding  layers, 
called  the  “outer  region”  (OR),  are  approximated  in  the  asymptotic  limit.  Approximate 
mode  equations  were  also  derived  by  matching  the  boundary  conditions  at  the  interface 
between  the  IR  and  the  OR,  whose  accuracy  can  be  improved  arbitrarily  by  including 
more  layers  in  the  IR.  For  the  sake  of  simplicity,  we  begin  our  analysis  of  air-core 
Bragg  fiber  by  matching  boundary  conditions  at  the  air  core  boundary  and  improve  the 
accuracy  of  our  results  later. 

Consider  the  generic  air-core  Braggg  fiber  given  in  Fig.  1,  the  mode  dispersion  rela¬ 
tions  depend  on  all  five  parameters:  n\,  n 2,  rcore,  r±  and  r2-  In  the  limit  of  large  radius, 


1In  Ref.  [10]  it  is  called  the  “core  region”,  but  we  decide  to  change  its  name  to  “inner  region”  so  as 
not  to  confuse  it  with  the  actual  fiber  core. 
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Fig.  1.  Schematic  of  an  air-core  Bragg  fiber.  The  fiber  cladding  consists  of  alternating 
layers  of  dielectric  media  with  high  and  low  refractive  indices. 

the  circular  Bragg  fiber  approaches  a  planar  Bragg  stack,  whose  band  diagrams  and 
mode  equations  can  be  found  in  standard  literature  [17].  Hence  for  low-loss  transmis¬ 
sion  it  is  necessary  to  have  the  Bragg  fiber  modes  to  fall  inside  the  TE  and/or  TM  band 
gap(s)  of  the  corresponding  Bragg  stack,  whichever  is  applicable. 

We  start  with  a  large  index  contrast  between  the  cladding  layers,  so  as  to  result  in 
a  large  band  gap.  We  choose  n\  =  4.6  and  ri2  =  1.5,  just  as  in  Ref.  [10]  and  we  try 
different  combinations  of  { rcore ,  r±,  r2j-  until  we  find  a  single-mode  frequency  window 
in  which  only  m  =  1  modes  with  negative  dispersion  exist,  where  m  is  the  azimuthal 
index.  Note  in  general  modes  in  a  cylindrically  symmetric  waveguide  are  designated  by 
two  mode  numbers,  one  radial  and  one  azimuthal.  In  this  paper  we  discuss  only  modes 
with  the  lowest  radial  mode  number  -  all  higher-order  modes  have  been  “pushed”  out 
of  the  band  gap  by  continually  shrinking  the  fiber  air  core. 

The  fiber  modes  can  also  be  classified  by  their  polarization  characteristics:  the  modes 
are  purely  TE  and  TM  for  m  =  0,  but  for  m  0  they  are  either  TE-like  or  TM-like 
(just  like  the  EH  and  HE  modes  in  conventional  fibers).  It  turns  out  that  the  TM  and 
the  TM-like  modes  with  the  lowest  radial  number  are  the  last  ones  to  remain  in  the 
band  gap  as  we  shrink  the  core  radius.  For  this  reason  only  these  modes  and  the  TM 
band  structure  (of  the  corresponding  Bragg  stack)  are  plotted  in  Fig.  2. 

In  Fig.  2  we  see  that  in  the  frequency  window  (0.705,  0.782)  only  the  m  =  1  mode 
is  present  and  form  a  band  in  the  fundamental  band  gap.  The  segment  of  the  band 
with  negative  curvature  or  dispersion  (i.e.  d2uj/d(5 2  <  0)  is  enclosed  in  the  shaded 
box  spanning  the  frequency  range  (0.762,  0.782).  Since  both  the  angular  frequency  u> 
and  the  propagation  wave  vector  /3  have  been  normalized  in  Fig.  2,  they  need  to  be 
“denormalized”  before  we  can  compare  results  with  previous  works.  For  instance,  if  we 
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Bragg  fiber  band  diagram 


Fig.  2.  Bragg  fiber  band  diagram  based  on  asymptotic  computations.  For  clarity 
only  the  TM  and  TM-like  modes  with  lowest  radial  number  are  plotted  here,  as  all 
other  modes  have  all  been  “pushed”  out  of  the  band  gap.  The  TM  band  structure  of 
the  corresponding  Bragg  stack  serves  as  the  background  and  the  dotted  line  shown 
is  the  light  line.  Boundary  conditions  are  matched  at  the  air  core  boundary.  The 
structural  parameters  are  given  as  follows:  n\  =  4.6,  n 2  =  1.5,  r\  =  10,  r 2  =  2,  and 
rCore  =  30,  all  in  normalized  units.  Note  in  the  frequency  window  (0.762,  0.782) 
the  m  =  1  band  not  only  exhibits  single-mode  behavior  but  also  possesses  negative 
dispersion  values. 

pick  the  normalized  operating  frequency  to  be  0.77(c/(ri  +  ^2)),  which,  when  converted 
to  1.55/rm  in  MKS  units,  gives  a  dispersion  value  D  of  roughly  -25,000  ps/(nm.km),  a 
25-fold  improvement  over  previous  results.  See  Fig.  3  for  a  plot  of  the  dispersion  value 
D  in  the  vicinity  of  the  operating  wavelength. 

No  optimization  technique  is  used  in  plotting  Fig.  2,  which  is  simply  the  result  of 
trial  and  error.  The  value  of  rcore  is  chosen  to  maximize  the  bandwidth  of  the  single¬ 
mode  window:  too  small  of  a  core  radius  would  push  even  the  m  =  1  mode  outside  the 
bandgap,  whereas  too  big  of  a  radius  would  introduce  too  many  higher-order  modes  into 
the  bandgap.  rq  and  r 2  are  also  chosen  based  on  similar  considerations.  As  a  general  rule, 
we  find  that  higher  r\jri  ratios  lead  to  higher  negative  dispersion  values.  This  however 
does  not  mean  that  we  should  let  rq/r2  increase  without  bound,  for  the  bandgap  shrinks 
to  zero  as  r\jr2  approaches  infinity  (Maximum  bandgap  occurs  when  the  quarter  wave 
condition  is  satisfied.  Any  deviation  from  it  results  in  a  reduced  bandgap).  This  trade-off 
between  single-mode  bandwidth  and  dispersion  value  D  gives  our  choice  of  rq  and  r2  in 
Fig.  2. 

3  Simulation  results 

We  have  already  used  the  asymptotic  algorithm  to  study  the  dispersion  properties  of  a 
Bragg  fiber  in  Sec.  2.  Now  we  use  finite  difference  time  domain  (FDTD)  simulations  to 
verify  the  validity  of  our  asymptotic  approach.  Note  our  inherently  3D  structure  can  be 
reduced  into  a  2D  one  by  extracting  out  the  exponential  3  dependence.  As  a  result,  if  we 
write  the  electric  field  as  E(x,y,t)exp(i/3z)  and  the  magnetic  field  as  H(x,y,t)exp(i/3z), 
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Dispersion  value  D  for  m=1  mode 


Fig.  3.  Dispersion  value  D  for  the  m  =  1  band  based  on  asymptotic  computations. 

a2  /o 

The  second  order  derivative  -??—£■  is  calculated  by  a  direct  differential  method  (i.e., 
no  curve  fitting  is  involved).  Fiber  structure  and  boundary  matching  surface  follow 
those  of  Fig.  2.  Shaded  area  (1526,  1566)  is  made  to  match  the  frequency  window 
(0.762,  0.782)  (also  in  Fig.  2)  by  choosing  the  normalized  operating  frequency  to  be 
0.77(c/(ri  +  7*2)),  which  translates  to  1.55/xm  in  MKS  units.  Note  at  A  =  1.55 /xra, 
dispersion  value  D  «  —25,000 ps/(nm.km). 

the  3D  Maxwell  equations  become: 

d 

(Vj.  +  i/3ez)  x  E(x,y,t )  =  -n0—H(x,y,t)  ,  (1) 

-  d  - 

(V_l  +i/3ez)  x  H(x,y,t)  =  e0e(x,y)  —  E(x,y,t)  ,  (2) 

where  is  defined  as  exj ^  +  ey-§^-  We  can  transform  these  2D  differential  equations 
into  a  set  of  finite  difference  equations  following  the  standard  finite  difference  time 
domain  procedure  [18].  At  the  boundaries  of  the  computational  domain,  we  use  the 
perfectly  matched  layer  boundary  condition  (PML)  [19,  20]  to  absorb  all  the  outgoing 
radiation.  For  the  2D  FDTD  equations  and  other  details  about  the  implementation  of 
the  algorithm,  the  reader  should  consult  Ref.  [21]. 

As  in  Sec.  2  we  consider  an  air-core  Bragg  fiber  with  the  following  parameters: 
v core  =  30,  Tii  —  4.6,  T\  =  10,  ri2  =  1.5  and  r2  =  2,  where  the  parameters  are  defined 
in  Fig.  1.  Twelve  cladding  pairs  are  used  in  our  simulations  to  reduce  the  radiation  loss 
to  approximately  0.2dB/km  [10]. 

To  verify  the  dispersion  curve  in  Fig.  2  with  FDTD  simulations,  we  first  give  the  fiber 
an  initial  field  and  let  it  evolve  in  time  in  accordance  with  the  Maxwell’s  equations,  with 
the  value  of  a  chosen  field  component  recorded  at  every  time  step.  We  then  proceed  to 
perform  a  fast  Fourier  transform  (FFT)  on  the  time  evolution  data,  and  each  confined 
mode  would  be  represented  by  a  peak  in  the  FFT  plot  in  frequency  spectrum. 

Both  the  asymptotic  and  the  FDTD  results  are  shown  in  Fig.  4,  where  the  asymptotic 
curve  is  simply  copied  over  from  Fig.  2.  We  find  a  good  agreement  between  the  two 
curves. 


#1376 -$15.00  US 
(C)  2002  OSA 


Received  June  18,  2002;  Revised  August  22,  2002 
26  August  2002  /  Vol.  10,  No.  17  /  OPTICS  EXPRESS  903 


Bragg  fiber  band  diagram  for  m=1  mode 


0.1  0.2  0.3  0.4  0.5  0.6 

P  (in  units  of  1/(r1+r2)) 


Fig.  4.  Bragg  fiber  band  diagram  for  the  m=l  mode.  The  asymptotic  curve  is  copied 
from  Fig.  2  and  the  FDTD  structure  is  defined  with  12  cladding  pairs. 

In  Fig.  5,  we  show  the  distribution  of  the  Hz  field  obtained  from  FDTD  calculations. 
The  frequency  and  propagation  constant  of  the  mode  are  respectively  uj  =  0.777 (c/ (r i  + 
r2))  and  f3  =  0.5(l/(ri  +  r2)).  Fig.  5  clearly  shows  that  the  guided  mode  has  a  radial 
number  l  =  0  (no  null  in  the  radial  direction  inside  the  core),  an  azimuthal  number 
m  =  1  and  a  significant  portion  of  the  field  concentrated  within  the  air  core. 

4  Observations  and  comments 

In  Sec.  2  we  found  a  frequency  window  in  which  only  the  m  =  1  mode  can  propagate. 
This  claim  was  made,  however,  based  on  asymptotic  calculations  that  matched  the 
boundary  conditions  at  the  air  core  boundary  (i.e.,  only  one  layer  in  the  IR).  Would  the 
claim  still  hold  true  as  we  include  more  layers  in  the  IR,  which,  as  Ref.  [10]  shows,  would 
improve  the  accuracy  of  the  solutions?  The  answer  is  an  unfortunate  ‘No’:  given  sufficient 
number  of  cladding  layers  we  have  found  numerous  higher-order  cladding  modes  in  the 
above-mentioned  single  mode  frequency  window.  One  such  higher-order  mode  is  plotted 
in  Fig.  6,  which  is  found  by  matching  the  boundary  conditions  at  the  outer  surface  of  the 
50th  layer  (i.e.,  the  IR  contains  50  layers).  This  m  =  4  band  cuts  through  the  frequency 
window  (0.762,0.782),  rendering  it  multi-mode.  The  m  =  1  curve  is  also  redrawn  under 
this  new  boundary  matching  condition  and  it  does  not  change  much  from  Fig.  2.  In 
fact  our  calculations  show  that  as  we  move  the  boundary-matching  surface  from  the 
air  core  boundary  to  the  50th  cladding  layer,  the  m  =  1  band  converges  very  quickly 
from  the  band  in  Fig.  2  to  the  one  in  Fig.  6.  Moving  the  boundary-matching  surface 
further  out  into  the  cladding  region  would  be  unnecessary  as  it  would  result  in  only  a 
minimal  correction  to  Fig.  6.  Thus  we  can  regard  the  m  =  1  band  in  Fig.  6  as  the  exact 
mode  solutions  and  work  out  the  “exact”  dispersion  value  D  at  the  same  frequency  of 
0.77(c/(ri  +  r2)),  which  turns  out  to  be  roughly  -20,000  ps/(nm.km). 

Just  as  a  side  note,  the  m  =  4  band  in  Fig.  6  is  not  the  same  as  the  one  in  Fig.  2. 
In  fact,  as  we  move  the  boundary-matching  surface  from  the  air  core  boundary  to  the 
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Fig.  5.  The  Hz  field  distribution  of  an  m  =  1  mode  based  on  FDTD  simulation.  The 
parameters  of  the  Bragg  fiber  are  given  in  the  caption  of  Fig.  2.  The  Bragg  cladding 
consists  of  12  cladding  pairs  and  the  whole  fiber  is  immersed  in  air.  The  frequency 
and  propagation  constant  of  the  mode  are  respectively  co  =  0.777(c/ {r\  +  r2))  and 
P  =  0.5(1/ (ri  +  r2). 


50th  cladding  layer,  the  m  =  4  mode  in  Fig.  2  gradually  moves  up  in  frequency  and 
eventually  disappears  from  the  bandgap.  In  other  words,  its  presence  in  Fig.  2  is  merely 
a  result  of  numerical  errors  inherent  in  our  asymptotic  calculations. 

Single  mode  transmission  would  seem  to  be  unrealizable  if  not  for  the  following 
considerations: 

(1)  Just  like  the  cladding  modes  in  conventional  fibers,  the  higher-order  cladding 
modes  in  Bragg  fibers  also  have  a  much  more  spread-out  field  profile  than  the  m  = 
1  modes,  which  translates  into  higher  loss  for  fibers  with  fixed  number  of  cladding 
layers.  This  is  also  manifested  in  the  fact  that  our  first  round  of  asymptotic  calculations 
(by  matching  boundary  conditions  at  the  air  core  boundary)  failed  to  detect  any  of 
the  higher-order  modes,  due  to  large  computation  errors  incurred  by  the  asymptotic 
approximations.  Theoretically  this  means  we  can  classify  the  Bragg  fiber  modes  into 
two  categories  based  on  their  confinement  characteristics:  those  of  truly  confined  air- 
core  modes  and  those  of  less  confined  cladding  modes.  In  Ref.  [10]  we  pointed  out  that 
Bragg  fibers  are  not  true  Bragg  structures  due  to  their  lack  of  translational  periodicity. 
Rather,  we  can  view  them  as  Bragg  structures  with  distributed  perturbations,  with  the 
amount  of  perturbation  decaying  to  zero  as  r  — >  oo.  Since  the  asymptotic  calculations 
assume  true  Bragg  confinement  after  the  boundary-matching  surface,  then  any  modes 
that  can  be  detected  by  matching  the  boundary  conditions  at  the  air  core  boundary  is 
confined  by  true  Bragg  reflections,  whereas  all  the  other  modes  require  at  least  some 
distributed  perturbations  to  achieve  confinement. 
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Bragg  fiber  band  diagram 
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Fig.  6.  Bragg  fiber  band  diagram  based  on  asymptotic  computations.  Just  as  in  Fig. 

2,  only  the  TM-like  bands  with  the  lowest  radial  number  are  plotted.  Boundary 
conditions  are  matched  at  the  50th  cladding  layer.  The  fiber  structural  parameters 
are  given  in  the  caption  of  Fig.  2.  Note  that  the  m  =  4  band  was  missing  in  Fig.  2. 

A  more  precise  comparison  is  made  in  Fig.  7,  in  which  the  field  profiles  of  both  the 
m  =  1  and  m  =  4  modes  are  plotted  at  the  operating  frequency.  We  find  that  the  m  =  1 
mode  is  much  more  confined  that  the  higher-order  mode. 

(2)  Modes  of  different  azimuthal  symmetries  are  orthogonal.  If  the  incoming  fiber  is 
properly  aligned  with  the  Bragg  DCF,  then  the  incident  m  =  1  mode  would  theoretically 
only  excite  the  m  =  1  mode  in  the  Bragg  fiber.  Of  course  in  practice  fiber  misalignment 
and  material  imperfections  would  never  allow  excitation  of  m  =  1  mode  alone  in  the 
DCF,  but  we  believe  it  could  still  be  a  contributing  factor  that  we  can  take  advantage 
of. 

In  our  FDTD  simulations  described  in  Sec.  3,  no  higher-order  cladding  modes  have 
ever  been  detected.  This  is  probably  due  to  both  large  radiation  losses  and  selective 
mode  excitation  by  the  initial  field,  confirming  the  validity  of  the  two  arguments  above 
(Given  enough  cladding  layers  and  a  proper  initial  field,  higher-order  modes  like  the 
m  =  4  mode  in  Fig.  6  can  be  found  using  FDTD  simulations  and  the  results  agree  well 
with  that  given  by  the  asymptotic  approach). 

We  want  to  point  out  again  that  this  paper  represents  only  a  theoretical  study  of 
the  dispersion  properties  of  Bragg  fibers  and  does  not  represent  a  design  of  a  realistic 
dispersion-compensating  fiber.  In  this  paper  we  only  aim  to  show  the  possibility  of 
achieving  large  negative  dispersion  through  Bragg  reflections.  The  choice  of  large  index 
contrast  in  the  cladding  layers  is  a  mere  convenience.  It  not  only  affords  us  a  large 
band  gap  to  work  with,  but  through  tight  mode  confinement  also  allows  us  to  verify 
our  asymptotic  calculations  within  a  reasonable  time  frame  using  FDTD  simulations.  In 
fact  our  calculations  show  that  the  large  negative  dispersion  value  obtained  is  a  direct 
consequence  of  Bragg  reflection  and  its  order  of  magnitude  does  not  depend  on  the  index 
contrast.  Also  note  the  material  dispersion  is  deliberately  left  out  of  our  calculations 
and  simulations,  as  it  depends  on  the  particular  choice  of  materials. 
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E  Field  Profile 


Fig.  7.  The  radial  profile  of  Eg  field  of  guided  Bragg  fiber  modes.  The  white  area 
is  the  fiber  core  and  the  shaded  layers  represent  the  cladding  layers.  All  structural 
parameters  are  given  in  the  caption  of  Fig.  2.  The  vertical  dashed  line  shows  the 
boundary  of  a  fiber  structure  with  12  cladding  pairs,  used  in  the  FDTD  simulations 
in  Sec.  3.  The  frequency  of  the  mode  is  to  =  0.77(l/(r*i  +7*2)),  the  chosen  normalized 
operating  frequency. 

5  Future  Work 

Can  we  build  a  realistic  dispersion-compensating  fiber  as  it  is?  The  answer  is  “No”. 
Though  Bragg  fibers  with  large  index  contrast  have  in  fact  been  made  [2],  with  the 
high-index  layer  made  of  Tellurium  (n  =  4.6),  such  a  high-index  material  suffers  from 
an  absorption  loss  that  is  generally  several  orders  of  magnitude  higher  than  that  of 
glass,  hence  makes  them  unusable  in  long-haul  communication  systems.  Besides,  the 
structural  parameters  of  the  Bragg  fiber  used  in  this  paper  are  not  compatible  with  the 
fiber  parameters  used  in  the  real  world.  For  example,  if  the  normalized  frequency  of 
0.77(c/(ri  +  7~2))  is  converted  to  the  operating  wavelength  A  =  1.55^tm,  then  the  fiber 
core  radius  rcore  works  out  to  be  roughly  0.47/im  (with  r\  =  0.16 /jm,  and  r2  =  0.032p,m 
respectively) ,  which  is  one  order  of  magnitude  smaller  than  that  of  a  conventional  fiber 
core,  hence  a  huge  splicing  loss  would  result.  One  possible  way  to  get  around  these 
problems  is  to  reduce  the  index  contrast  between  the  cladding  layers.  The  problem, 
however,  is  that  Bragg  fibers  of  low  index  contrast  require  a  large  number  of  cladding 
layers  for  mode  confinement.  For  example,  Bragg  fibers  with  an  index  contrast  of  0.1  in 
the  cladding  layers  require  about  1,000  cladding  pairs  to  achieve  0.2  dB/km  loss  [10], 
which  poses  a  serious  problem  to  fiber  manufacturers.  Research  on  reducing  the  number 
of  cladding  layers  required  for  low-index  fibers  would  be  the  subject  of  a  future  paper. 

As  a  final  remark,  we  note  that  in  this  paper  we  have  only  concerned  ourselves  with 
the  dispersion  value  D,  and  ignored  the  other  important  figure  of  merit  for  dispersion 
compensation,  namely,  the  dispersion  slope.  This  is  so  because  (1)  we  believe  that  in  a 
realistic  Bragg  DCF  with  low  index  contrast,  the  single-mode  bandwidth  would  be  so 
small  (due  to  a  very  narrow  bandgap)  that  multi-channel  compensation  is  unlikely.  (2) 
there  is  actually  a  frequency  at  which  the  dispersion  D  “bottoms  out”,  i.e. ,  dispersion 
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slope  is  zero  at  this  point.  For  the  Bragg  fiber  in  Fig.  2,  this  zero-slope  point  is  located 
well  above  the  the  single-mode  window,  and  hence  is  not  shown  in  Fig.  3.  Luckily,  we 
have  observed  a  sizable  downward  shift  in  its  frequency  as  we  reduce  the  index  contrast, 
which  should  work  in  our  favour  when  we  develop  the  DCF  with  low  index  contrast. 
Once  again,  the  details  will  be  published  in  a  future  paper. 

6  Conclusions 

In  this  paper,  we  studied  the  dispersion  property  of  the  m  =  1  mode  in  an  air  core  Bragg 
fiber,  which  lends  itself  well  to  dispersion  compensation  owing  to  reduced  absorption 
loss  and  nonlinear  effects.  Based  on  asymptotic  calculations  we  showed  large  negative 
dispersion  values  can  be  achieved  and  the  accuracy  of  the  algorithm  has  been  verified 
by  FDTD  simulations. 
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